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Section  1 

INTRODUCTION  AND  BACKGROUND 


The  Acoustic  Bottom  Interaction  contract 
( N00014-79-C-0370)  consists  of  three  tasks,  which  are 
summarized  as  follows: 

1.  Develop  a  geophysical  bottom  model  using  NADC 
bottom-loss  data  and  other  available 
acoustical  and  geophysical  data,  from  which 
low  frequency  bottom-loss  versus  grazing  angle 
and  frequency  can  be  derived. 

2.  Examine  the  NADC  thin-sed iment  bottom-loss 
data  to  ascertain  the  dominant  loss 
mechanisms. 

3.  Investigate  the  feasibility  of  identifying  the 
limits  on  the  likely  range  of  spatial 
variability  of  small-scale  structure  in 
unconsolidated  sediments,  and  determine  the 
feasibility  of  using  existing  random-medium 
propagation  models  to  relate  coherence  of  the 
acoustic  field  to  this  spatial  variability. 
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Section  2 

GEOPHYSICAL  BOTTOM  MODEL 


SAI  has  developed  a  geophysical  bottom  model  based 
on  NADC  bottom  loss  data  and  available  geophysical  data. 
The  word  model  is  used  in  two  senses.  In  the  first 
instance,  the  parameters  necessary  to  specify  sediment 
acoustic  properties  are  identified.  Following  that, 
specific  models  in  the  form  of  a  numerical  data  base 
are  evaluated  and  extrapolated  over  regions  covering  the 
northern  hemisphere.  The  result  is  called  a  geophysical 
model  since  the  acoustical  volume  properties  of  the  sedi¬ 
ment,  such  as  velocity  and  attenuation  profiles,  are  used 
and  are  rela.ted  to  the  geological  history  of  the  sediment. 

A  set  of  ten  geophysical  parameters  are  identified 
that  give  good  agreement  between  predictions  and  NADC  bottom 
loss  data.  This  is  accomplished  using  a  ray  theory  model 
which  includes  the  effects  of  discontinuities  and  sound 
speed  and  attenuation  profiles.  For  purposes  of  comparing 
with  data,  the  known  processing  artifacts  due  to  Lloyd's 
mirror  and  octave  processing  are  also  modeled.  The  ten 
parameters  are:  sound  speed  ratio  of  the  water/sediment 
interface,  two  parameters  describing  the  sound  speed  pro¬ 
file,  two  parameters  for  the  attenuation  profile,  bulk 
sediment  density,  thin  surface  reflection  with  an  associated 
density  and  thickness,  and  a  sediment  thickness  and  basement 
reflection  coefficient.  A  brief  discussion  of  each  para¬ 
meter  follows. 
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The  sediment  sound  speed  profile  is  based  on 
Hamilton's  (1979)  average  estimates  of  profiles  for  sedi¬ 
ments  that  have  terrigenous,  calcareous,  or  silicious 
content.  These  profiles  are  modified  somewhat  by  the 

[  fact  that  sound  speed  varies  at  the  sediment  top,  and  is 

governed  by  a  sound  speed  ratio,  i.e.,  the  ratio  of 

\  sedimc nt-to-water  sound  speeds  at  their  interface.  These 

*  - 

ratios  are  estimated  from  NADC  data.  Our  estimates  of 
the  ratios  agree  fairly  well  with  estimates  based  on  our 

best  understanding  of  local  geology  and  Hamilton's  (1980) 
canonical  estimates. 

The  sediment  attenuation  profile  increases 
linearly  with  depth.  It  is  the  least  well  understood 
acoustic  property  of  sediments,  but  this  form  of  profile 
is  consistently  required  to  fit  NADC  bottom  loss  measure¬ 
ments  over  the  full  range  of  geological  and  physiographic 
environments. 

Sediment  density  and  thin  layer  density  and 
thickness  are  an  approximation  of  the  actual  reflection 
process  in  the  sediment.  In  reality,  the  sediment  reflec¬ 
tion  process  cannot  be  modeled  as  a  two  fluid  problem. 
Sediments  are  generally  layered  due  to  different  epochs  of 
geological  deposition.  This  is  particularly  true  in  abyssal 
plains.  They  have  many  small  acoustic  reflections  near 
their  surface  that  contribute  to  the  reflected  signal,  and 
can  produce  strong  reflections,  which  are  extended  in  time. 
When  the  signal  is  processed  using  total  energy,  multiple 
layers  can  be  lumped  into  a  single  layer  to  give  a  good 
approximation  of  the  reflection  process.  This  is  verified 
independently  by  workers  at  ARL  in  Austin,  Texas.  Sedi¬ 
ment  thicknesses  are  estimated  on  a  world  wide  basis  and 
presented  in  the  form  of  contour  maps.  This  data  is 
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interpolated  to  produce  sediment  thicknesses  in  a  sixth 
of  a  degree  data  base.  These  results  are  based  on  partial 
results  from  Lamont  Doherty  Geological  Observatory  and 
bathymetry  information.  Basement  reflectivity  is  modeled 
with  a  single  reflection  coefficient,  and  is  sufficient  for 
thick  seaiment  areas. 

Geological  information  is  used  to  support  the 
extrapolation  of  results  of  the  acoustical  measurements 
at  specific  sites  to  the  full  ocean  bottom.  This  is 
derived  from  a  variety  of  sources,  predominantly  from 
reports  of  the  Deep  Sea  Drilling  Project  (DSDP).  Other 
information  of  a  general  nature  is  available  in  the  standard 
geological  literature  in  well  explored  areas  such  as  the 
Bermuda  Rise  and  the  Greenland  UK  gap.  Information  relating 
to  the  Navy  long-range  propagation  problem  is  organized  by 
bottom-loss  region  for  the  Northern  hemisphere  south  of 
70°  north.  This  consists  mainly  of  grain  size  and  chemi¬ 
cal  composition  of  bulk  and  surface  sediments  for  each  area 
as  well  as  physiographic  province  and  history  of  deposi¬ 
tion.  Grain  size  and  chemical  composition  is  used  to 
estimate  properties  of  sound  speed  iatio  and  depth  depend¬ 
ence  of  the  sound  speed  profile.  This  together  with  the 
history  of  sediment  deposition  is  used  to  extrapolate 
measured  acoustic  properties  to  areas  where  no  direct 
acoustic  measurements  are  available. 

The  geological  analysis  has  been  a  very  large 
part  of  this  task  since  geological  information  is  exten¬ 
sive  and  scattered  throughout  the  literature.  Further¬ 
more,  where  it  is  organized  it  is  generally  to  support  a 
geological  hypothesis  such  as  sea  floor  spreading.  In 


this  case,  acoustic  information  is  not  highlighted  and  in 
some  cases  is  suppressed.  It  should  be  emphasized  here  that 
while  the  bulk  of  geological  exploration  has  dealt  with  the 
average  properties  of  deep  sediments,  Navy  acoustic  propaga¬ 
tion  is  governed  by  the  detailed  properties  of  near  surface 
sediments.  The  current  project  provides  the  first  world 
wide  compilation  of  geological  information  relevant  to  the 
Navy  propagation  problem.  Together  the  acoustical  measure¬ 
ments  and  geological  information  have  led  to  the  identifica¬ 
tion  of  the  relevant  acoustic  parameters,  and  their  correla¬ 
tion  with  geological  properties  of  sediments.  This  has  led 
to  the  ability  to  extend  our  estimates  of  acoustic  para¬ 
meters  to  areas  where  geological  information  is  known.  On 
this  basis,  a  consistent  set  of  geoacoustic  model  parameters 
have  been  assigned  for  deep  ocean  areas  of  the  northern 
hemisphere. 
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Section  3 

THIN  SEDIMENT  LOSS  MECHNANISMS 


Task  2  is  comprised  of  two  sub  tasks  that 
examine  NADC  thin-sediment  bottom-loss  data  in  the  context 
of  a  mathematical/computer  model.  The  model  attempts  to 
reconcile  the  data  with  known  sediment  refraction  and  loss 
mechanisms  and  a  model  for  rough  basement  random  scattering. 
The  data  are  used  to  initially  estimate  geophysical  and 
model  parameter  inputs  for  the  specific  ocean  region  and 
measurement  configuration  where  the  data  were  obtained. 
These  data  characteristically  show  modest  to  high  bottom- 
losses  at  nearly  all  angles  and  frequencies,  with  loss 
increasing  with  frequency  but  nearly  independent  of  grazing 
angle.  Any  proposed  model  must  conform  to  these  data 
characteristics  in  order  to  correctly  ascertain  the  dominant 
loss  mechanisms. 

The  first  sub  task  concerns  a  rough  surface 
scattering  model  that  would  be  appropriate  for  basement 
scattering  and  consistent  with  the  data  behavior.  A 
detailed  description  of  such  a  mathematical  model  follows 
below.  The  second  sub  task  concerns  the  computer  imple¬ 
mentation  of  the  rough  surface  basement  scattering  model 
integrated  with  an  appropriate  water/sediment  refraction  and 
loss  model.  A  detailed  description  of  this  implementation 
follows  the  mathematical  model  description. 
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3.1 


SOUND  SCATTERING  FROM  A  ROUGH  OCEAN  BASEMENT 


To  analyze  the  effect  of  a  rough  ocean  basement  on 
sound  propagation  in  the  ocean,  the  following  model  is 
used . 

The  basement  is  described  as  a  random  rough 
surface,  with  a  zero  mean  Gaussian  distribution  of  heights. 
Further,  for  simplicity,  the  surface  is  assumed  wide  sense 
stationary,  i.e.,  the  correlation  function  W  between  two 
points  of  the  surface  depends  only  on  the  distance  between 
the  points;  it  is  also  isotropic,  i.e.,  independent  of 
direction.  The  surface  is  then  characterized  by  its  root- 
mean-square  (rms)  height  a,  its  correlation  distance  i , 
and  its  power  spectrum  S(k),  the  Fourier  transform  of  the 
correlation  function,  the  rms  slope,  -yA  a/l  ,  of  the  surface 
is  then  given  by  y  =  v -W"(0) .  The  distribution  of  surface 
slopes  is  again  Gaussian,  with  zero  mean. 

The  assumption  of  a  Gaussian  surface  allows  the 
use  of  a  composite  surface  theory  due  to  Brown  (1980).  The 
surface  is  decomposed  into  two  (Gaussian)  surfaces,  one 
containing  large  scale  features  (small  wave  number),  and 
superimposed  on  that,  the  second,  describing  small  scale 
(ripple)  effects. 


The  purpose  of  this  decomposition  is  as  follows. 
Scattering  from  large  scale  surface  features  is  dominated  by 
"glints",  i.e.,  specular  reflection  of  the  incident  plane 
wave  from  properly  oriented  facets  of  the  large  surface. 
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Contributions  to  the  average  received  intensity  of  this 
effect  is  obtained  using  the  tangent  plane  (Kirchhoff) 
method,  (geometrical  optics  solution).  Small  surface 
features  produce  a  diffuse  scattering,  and  the  intensity 
contribution  is  obtained  as  a  first  order  perturbation  of 
the  geometrical  optics  solution.  Since  both  surface 
features  are  independent,  the  total  average  intensity  is  the 
incoherent  sum  of  these  two  contributions,  as  shown  by 
Brown  (1980). 


The  division  between  large  and  small  scale  sur¬ 
faces  is  in  terms  of  a  surface  wave  number  Kj).  The  large 
surface  has  a  spectrum  Sl  (K)  defined  by 


SL  (K) 


S(K),  0<K<KDj 
0,  K  >  Kd 


while  the  small  scale  surface  spectrum  Ss(K)  is  given  by 


Ss  (K)  = 


o.  0  <  K  >  Kt 


,S  (K),  K>Kd 


Kd  is  determined  as  follows.  The  rms  height  os  of  small 
surface  features  is  defined  by 


a  =  2tt 
s 


f 

J  Kt 


k  S(k )dk 


where  Kd  is  chosen  so  that  the  Rayleigh  number  of  small 
surface  features  is  small,  i.e.,  4k02as^<<l,  and  k0  is  the 
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wave  number  of  the  incident  plane  wave.  This  requirement 
determines  KD.  (In  practice,  we  set  4kQ2as2  =  .13,  and 
solved  for  Kp). 


Having  determined  Kq,  the  rms  slope  Y  l  of  large 
surface  features  is  given  by 

Kr 


y L  =  -W"(0)  =  2 r 


J 


LD 


k  S(k)dk 


The  Gaussian  slope  distribution  of  large  surface  features  is 

2 

given  in  terms  of  Y l  . 


The  mean  scattered  intensity  from  large  surface 
features  is  given  by  (see  Bass,  articles  20  and  24) 


|c  I2 


■/- 


i  v  L2_al-. 

4Ri2  R2  4 


dS 


the  overbar  indicates  "average' 


is  the  ensonified 


area  of  the  ocean  basement,  and  if  K^,  K2  are  the  incoming 
and  outgoing  wave  vectors  at  the  basement,  then  "q~ =  -  Kg. 

,  Kg  have  norm  k0 ,  and  incident  and  azimuthal  angles 
>  <f>o >  and  b  ,  <f> ,  respectively.  The  azimuthal  angles 
are  measured  relative  to  the  plane  which  contains  the 
source  and  receiver,  and  is  normal  to  the  ocean  floor. 


q2  =  1*1  =  2k 


1  +  cose  COS0O  -  Sine  sine0  cos  ( <P0  - 


(K1  “  K2)z  =  kQ(cos  6  +  cos  e  ). 


']■ 
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y  =  -qi _/qz  defines' the  angle  of  specular  reflection  from 
source  to  receiver,  where  q =  (qx,  qy).  R1  and  R2  are 
distances  along  the  ray  paths  from  source  and  receiver 
respectively,  to  the  basement.  V  is  the  basement  reflection 
coefficient  (set  equal  to  one  in  this  preliminary  study). 


we  is  the  effective  distribution  of  surface 
slopes  of  large  scale  surface  features  in  the  presence  of 
shadowing.  Following  Bass  (articles  23  and  24),  and 
Sancer  (1969) 


w 


|y  -  -qL^=  w  /y  =  R  (0,  0Q)  , 


where  w(  Yx  Yy)  is  the  slope  density  function  of  large 
surface  features  given  by 


w  (y  ,  y  )  = 

x  ’  y 


exp 


•(y  2  +  y  2) 
v  1  x  1  y  ' 


27ryT“  \  2  Y] 

R(9,0O)  is  a  correction  due  to  shadowing,  given  by 


R(0,0O)  = 


1  +  C^d)  +  C ( 6 


•>] 


-1 


/2yf \  !/2  /~cot2e 1 

with  2C(0  *  )  =  I  tan0  '  •  exp  ( - 5 - 


2y 


L  ' 


-  erfc 


cot  0  1 

m t)1'2 


(Strictly  speaking,  the  above  formula  for  R  is 
correct  if  $  0-  4>  (the  azimuthal  difference)  is  near  tt  .  For 
back  scattering,  (  4>0-  <f>~0),  the  term  C(  0 )  +  C(  0O)  should 


be  replaced  by  C(max  (  8 ,  6 0 ) ) .  In  the  absence  of  any 
shadowing  theory  providing  a  transition  between  <j>  0- <P  ~  tt  and 
<f>0- <l>=:0 ,  it  was  decided  to  use  the  abo^'e  formula  for  all 
The  difference  is  expected  to  be  slight.) 

The  perturbation  term  due  to  diffuse  reflection 
from  small  surface  features  is  of  the  form  (see  Bass, 
articles  8  and  24,  and  Brown), 

qcT:  T)  R(A.  eD)  r(q,  -,L>  Ss)  ds 

o 

where  S0,  R* ,  R2»  kQ,  and  R($,<J>0)  are  as  before. 

=  koT:  TtJ  =  ko"T,  and  .  T )  is  given  b>  •  >'<»»•  \v\  =  1), 

q  k  .  6)  =  (1-  dx.  g^)2  =  {^l-sin6sin6o  cos  (0o-i>)J^ 

r(q,  yl,  Ss  )  is  a  convolution  integral  involving  the 
small  spectrum  Ss  and  the  slope  density  function  for  large 
surface  features  with  limits  involving  the  effects  of 
shadowing . 


For  this  preliminary  study,  F ( q ,  y  l  >  Ss)  is 
approximated  by  2  tt  Ss(  |q|  ).  (Brown  provides  a  rationale, 
using  results  of  Stogryn,  for  including  the  shadowing 
function  R  within  the  perturbation  integral.  The  physical 
notion  is  chat  diffuse  scattering  cannot  occur  if  the  small 
scale  surface  features  are  shadowed  by  large  scale  features. 
This  approximation  for  T  may  not  be  adequate  if  Y  l2  qz2 
becomes  too  large. 
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then 


The  total  average  intensity  at  the  receiver  is 


!  u  1 2+  I  u  I  2  =  f  K(6,e0>4>0-<{,la,Y,S,k0)dS, 

S  R1R2 
o 

where  the  kernel  K  has  been  implicitly  given  above. 


The  features  of 


1  a  1  +  a  w  ^  n  A  /N  4-  V»  rs  ^  ^  U  A  + 

kj  j.  u  n  jl  o  i/O  oui  xavc  ujcui  y  tiia  v 


makes  this  model  appealing  are  as  follows: 


First,  as  shown  by  Brown,  in  the  absence  of  large 

scale  features,  one  recovers  the  classical  solution  (see 

Rice)  for  scattering  from  a  slightly  perturbed  planar 

2 

surface.  In  the  limit,  asY  -*~0,  the  integral  produces 
the  specular  reflection  correctly.  Numerically,  the 
approximation  of  the  integral  over  SQ  is  accurate,  using  a 
2000  foot  mesh  size,  down  to  a  mean  slope  of  about  8°. 
Below  that,  a  smaller  mesh  size  over  SQ  .is  needed  to 
accurately  evaluate  the  specular  contribution. 

2 

Second,  the  choice  of  Kq>  and  therefore  Y^  , 
depends  on  K0.  This  introduces  a  frequency  dependence  in 
the  geometrical  optics  solution,  an  effect  which  has  been 
observed  in  radar  studies  of  the  moon  (see  Evans  and 
Hagfors).  Of  course,  the  usual  geometrical  optics  solution 
is  a  high  frequency  limit,  and  is  necessarily  independent  of 
k0 . 


Third,  following  Lynch  and  Wagner  (1970a,  1970b) 
the  kernel  has  been  checked  for  energy  conservation.  This 
requires  that  for  all  8  ,  0  Q,  o,  Y,S,  kQ, 


r  \  r2v 

J  I  K(e,eo.  $o-4,0,Y,Sfko)  Sine  d<f>d0  =  cosr0. 
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A  computation  has  been  made,  and  the  result  is  correct, 
up  to  a  factor  near  1,  at  most  1.18.  While  negligible,  the 
kernel  has  been  re-normalized  by  a  factor  of  (1.15)~1. 

It  should  be  noted  that  the  above  shadowing  theory 
is  a  single  scatter  theory  only,  containing  no  correction 
for  multiple  scatters.  The  present  theory  is  adequate  in 
this  respect  except  for  very  rough  surfaces  (mean  slopes 
above  40°),  and  incident  angles  near  zero. 

3.2  COMPUTER  IMPLEMENTATION 


...... 


The  rough  basement  scattering  integral  for  total 
average  intensity  given  above  is  one  part  of  a  complete 
environmental  acoustic  model  implemented  on  a  computer.  A 
complete  model  requires  that  the  incident  intensity  at  the 
sediment/basement  interface  and  at  the  receiver  correctly 
account  for  spreading  loss  and  refraction  in  the  water 
column  and  refraction  and  attenuation  in  the  sediment. 
Since  source-to-receiver  distances  are  not  great,  water 
column  attenuation  is  not  significant.  Also,  the  incre¬ 
mental  effects  of  travel  time  and  spreading  loss  due 
to  the  portion  of  the  ray  path  within  the  sediment  are 
small  enough  to  be  ignored. 


Spreading  loss  and  refraction  in  the  water  column 
(to  the  water/sediment  interface)  are  computed  using  a  ray 
trace  program  with  the  same  source/receiver  depths  and  sound 
speed  profile  used  by  NADC  and  measured  during  the  actual 
sea  tests.  Best  estimates  of  sediment  properties  are  ob¬ 
tained  from  historical  records  and  other  analyses,  however 
sediment  thickness  is  estimated  from  NADC  bottom-loss  data 


(the  procedure  used  is  discussed  in  detail  below),  and  a 
constant  sound-speed  gradient  is  also  estimated.  The 
basement  scattering  kernel  has  been  described  in  detail 
above.  Integration  over  the  scattering  region  is  approxi¬ 
mated  by  a  summation  over  a  grid  representing  the  ensoni- 
f ied/scattering  region. 

The  first  step  in  the  simulation  ic  to  compute  the 
e-nsonif ied/scattering  region  at  the  water/sediment  interface 
over  a  range  of  grazing  angles  from  zero  to  90  degrees.  A 
family  of  rays  are  traced  from  the  ocean  bottom  to  source 
and  receiver  depths,  allowing  up  to  one  surface  bounce.  By 
reciprocity  these  would  be  the  same  rays  from  source/ 
receiver  depths  to  the  bottom.  A  table  is  generated  that 
records  for  each  ray  traced:  bottom  grazing  angle,  spreading 
loss,  travel  time,  arrival  angles,  and  horizontal  range  from 
source/receiver  to  the  point  of  bottom  contact.  A  smooth 
pressure  release  surface  is  assumed. 

These  tables  provide  the  basic  input  needed  for 
generating  the  acoustic  field  incident  on  the  bottom  and 
scattered  field  to  the  receiver  for  as  many  source-to- 
receiver  distances  as  required  to  compute  bottom-loss  as  a 
function  of  grazing  angle.  The  total  ensonified  area  is 
given  by  a  circle,  centered  at  the  source  projection  on  the 
bottom,  of  radius  equal  to  the  horizontal  range  to  zero 
grazing  angle  for  a  single  surface  reflected  ray  leaving  the 
source  in  an  upward  direction.  Similarly,  the  scattering 
region  to  the  receiver  is  computed  as  a  circle  centered  at 
the  receiver  projection  on  the  bottom.  For  a  given  source/ 
receiver  distance  these  two  circles  intersect  if  the  dis¬ 
tance  is  less  than  the  sum  of  the  two  circle's  radii,  their 
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intersection  being  the  bottom  region  from  which  scattered 
incident  energy  can  reach  the  receiver.  A  uniform  grid  is 
then  set  up  to  encompass  this  region  and  all  tabulated  ray 
trace  values  are  interpolated  to  these  grid  points.  This 
new  table  now  serves  as  a  basis  for  computing  basement 
scattering  for  each  source/receiver  geometry. 

Before  this  is  done,  however,  a  sediment  is 
introduced  as  follows.  For  the  NADC  test  site  used,  a 
water /sediment  sound  speed  ratio  and  sediment  sound  speed 
gradient  are  estimated.  By  examining  NADC  bottom-loss  data, 
one  may  also  deduce  a  lower  grazing  angle  cut  off  below 
which  bottom  loss  rapidly  drops  off.  Assuming  rays  incident 
below  this  cut  off  angle  do  not  reach  the  basement  and 
therefore  do  not  reach  the  receiver,  a  sediment  thickness 
can  be  estimated.  Again,  using  NADC  bottom-loss  data  at 
large  grazing  angles  a  frequency  dependent  loss  is  estimated 
using 

Loss  =  BRL  +  2  a  W 

where  BRL  =  basement  reflection  loss  constant 

a  =  attenuation  coefficient/unit  length/frequency 
l  =  one  way  sediment  ray  path  length 
f  =  frequency 

Since  &  can  be  computed  and  f  is  known,  the  data  are  used 
to  solve  for  BRL  and  a ,  which  is  then  used  to  calculate 
sediment  loss  for  each  refracting  ray. 
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The  equation  for  total  average  received  intensity, 
given  in  Section  3.1  is  modified  to  account  for  an  inhomo¬ 
geneous  ocean  by  replacing  (1  R  i 2  R  2  2  )  by  ( I  ^  1 2 )  , 

(which  have  been  tabulated),  at  the  basement,  and  summing 
over  the  appropriate  bottom  grid  „jing 


V' 

Lj 


ilA2 


,<j>G  -  <b.o  ,Y  ,S  ,kG)  A  A 


subject  to  any  angle  or  travel  time  constraints  that  may  be 
appropriate.  At  each  grid  point  the  effects  of  the  four 
dominant  bottom  bounce  ray  paths  are  computed,  i.e.,  source, 
bottom  bounce,  receiver;  source,  surface,  bottom  bounce, 
receiver;  etc. 

A  time  history  of  arrival  intensities  is  also 
computed  using  the  specular  path  travel  time  (which  is 
minimum  bottom  bounce  path  travel  time)  as  reference  and  a 
specified  time  increment  bin.  Hence,  rays  arriving  within 
each  time  increment  are  subsummed  and  tabulated  along  with 
such  other  useful  related  information  as  average  grazing 
angle,  average  azimuth  angle,  number  of  rays  arriving  within 
the  interval,  etc.  These  data  may  also  be  computer  plotted, 
if  desired,  to  present  visual  synopses,  which  are  useful  for 
understanding  and  relating  scattering  phenomena  to  the 
various  geophysical  mechanisms,  e.g.,  roughness  parameters, 
sediment  parameters,  etc.  This  type  of  output  allows  a 
quasi-quantitative  comparision  with  available  time  history 
data. 
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A  limited  amount  of  data  have  been  used  for 
program  check  out,  and  several  rough  surface  spectra  have 
been  tried  for  the  purpose  of  determining  qualitative 
agreement  with  data  and  the  degree  of  surface  roughness  and 
sediment  properties  required. 

Work  is  continuing  in  this  area  under  a  follow-on 
contract  entitled  "Thin  Sediment  Bottom  Interaction  Measure¬ 
ments".  The  objective  of  this  contract  is  to  apply  the 
computer  code  with  appropriate  bottom  acoustic  properties 
and  scattering  kernel  to  interpret  a  set  of  acoustical 
data. 


fcf*  j  ' 
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Section  4 

THICK  SEDIMENT  COHERENCE 

Two  issues,  which  have  been  identified  as  having 
high  priority  to  the  Bottom  Interaction  Programs,  (Hanna  and 
Hawker),  (Martin),  are  the  rate  of  lateral  variation  of 
sediment  properties,  in  different  oceanic  regions  and  the 
degree  of  coherence  of  bottom-interacting  energy  that 
propagates  through  such  sediments  and  is  scattered  by 
sediment  inhomogeneities.  Of  fundamental  interest  is 
whether  this  bottom-interacting  energy  is  useful  for  signal 
processing  systems. 

The  final  report  for  this  task  (Dozier,  1981) 
(cf.  Enclosure)  shows  the  feasibility  of  using  an  extension 
of  an  existing  random-medium  propagation  model  (Dozier  and 
Tappert,  1978)  to  relate  coherence  of  the  acoustic  field  to 
lateral  variability  of  sediments.  The  approach  is  a 
coupled-mode  theory  valid  in  a  low-frequency  limit;  of 
course,  low  frequencies  are  precisely  where  acoustic  inter¬ 
action  with  sediment  is  important.  A  significant  complica¬ 
tion  not  included  in  the  earlier  work,  (ibid.)  however,  is 
the  inclusion  of  a  sediment  model,  with  density  and  attenua¬ 
tion  as  well  as  sound  speed.  Rutherford  and  Hawker  (1978) 
have  shown  that  sediment  density  gradients  play  a  very  minor 
role  in  bottom  reflection  loss;  hence,  we  choose  a  layered 
sediment  model  in  which  density  is  constant  within  each 
layer.  Sediment  sound  speed  is  then  modeled  as  varying 
randomly  about  a  mean,  deterministic  profile  in  depth. 
Although  in  fact,  unlike  the  earlier  case  of  internal  waves 
in  the  water  column  (Dozier,  JASA  63,  1978),  the  sediment 
sound  speed  is  deterministic,  we  lack  detailed  subbottom 
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data  and,  even  if  it  were  available,  it  would  be  too  complex 
to  be  useful.  For  our  purposes  the  scale  length  of  hori¬ 
zontal  variation  is  the  important  feature,  and  is  the  key 
ingredient  of  our  random  model. 

Thus,  we  derive  a  random  coupled  mode  equation 
which  describes  statistically  the  coupling  between  the 
deterministic  normal  modes  of  our  mean  sound  speed  profile, 
given  attenuation  and  piecewise  constant  density.  (Lateral 
variations  in  density  and  attenuation  can  be  included  also, 
but  are  relatively  unimportant  and  unnecessarily  complicate 
the  coupled  mode  equation.)  Unlike  the  previous  reference 
however,  the  coupling  coefficients  are  now  complex  (because 
of  attenuation),  and  are  neither  symmetric  nor  Karmitian. 
The  coupled  mode  equation,  therefore,  losses  its  energy- 
conserving  property  as  expected,  due  to  attenuation. 

The  next  step  is  averaging  the  coupled  mode 
equations  to  obtain  a  master  equation  for  the  mean  acoustic 
intensity.  Although  non-conservation  is  a  significant 
complication,  averaging  can  still  be  done.  The  master 
equation's  coefficient  matrix  is  no  longer  symmetric, 
and  modal  energies  do  not  approach  equipartition . 

Next,  the  averaging  is  extended  to  include  a 
two-point  statistic  of  horizontal  coherence  transverse  to 
the  direction  of  propagation  and  for  separation  distances 
much  less  than  the  source  range.  This  is  done  in  a  quasi¬ 
static  manner  by  treating  horizontal  distance  as  a  para¬ 
meter,  just  as  time  was  treated  earlier,  (Tappert  and 
Dozier,  1978),  and  results  is  a  master  equation  for 
horizontal  coherence. 


A  question  now  is,  when  iz  the  averaging  mathe¬ 
matically  valid?  It  is  done  using  a  perturbation  technique, 
but  without  precise  knowledge  of  the  horizontal  distance 
parameter  or  of  the  domain  of  validity.  The  distance 
parameter  is  essentially  determined  by  the  product  of  rms 
sound  speed  variation  and  correlation  length  of  the  varia¬ 
tion  (Dozier,  JASA  63,  1978).  Typical  rms  variations  found 
by  Matthews  (1980)  are  10%  (relative  to  the  mean)  for 
calcareous  sediments  and  20%  for  tu-bidite  sediments. 
Houtz  (1980)  finds  (experimentally)  standard  deviations 
mostly  from  10%  to  15%  over  large  ?rsas  of  the  North 
Atlantic.  These  relative  variations  are  roughly  10  times 
the  relatiave  variation  of  sound  speed  due  to  internal 
waves.  On  the  other  hand,  Tyce,  et.  al . ,  (1980),  have  found 
high  lateral  variability,  suggesting  that  the  correlation 
length  of  the  random  sound  speed  function  is  much  less  than 
the  11  km  distance  for  internal  waves  in  Dozier  (1981). 
Since,  as  noted  above,  it  is  the  product  of  correlation 
length  and  rms  sound  speed  variation  that  essentially 
determines  the  perturbation  parameter,  we  expect  a  domain  of 
validity  comparable  to  that  found  for  internal  waves,  i.  e., 
up  to  perhaps  400  Hz.  This  includes  most  of  the  low- 
frequency  regime  where  bottom  interactions  are  important. 

Until  more  bottom  data  becomes  available,  the 
primary  value  of  this  model  is  to  bound  the  coherence 
problem.  For  lateral  variability  within  the  model's  esti¬ 
mated  domain  of  validity,  a  numerical  implementation  will 
yield  the  associated  loss  of  cohere  ce.  Comparing  model 
results  for  different  scales  of  lateral  variability  will 
illustrate  trends  (e.g.,  loss  of  coherence  with  increasing 


rms  sound  speed  variation,  with  correlation  length  held 
fixed;  and  vice  versa)  and  point  out  worst  cases.  Cer¬ 
tainly,  for  many  real  ocean  environments  the  model  should 
predict  whether  a  significant  loss  of  coherence  is  likely. 

Unfortunately,  numerical  implementation  of  this 
model  was  beyond  the  scope  of  this  contract,  but  we  hope  to 
be  able  to  do  this  in  future  BIP  work. 
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